Analysis of the conduction heat transfer in cantilevers under steady state cryogenic 

conditions 
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An accurate analysis of the conduction heat transfer in a cryogenic flask is made and some useful 
formulae are derived. Taking into account the temperature dependence of conductivity and tensile 
strength of the supporting rods for a helium cryostat, these formulae may provide more exact results 
than the the formulae based on simpler models. This allows the design of the supporting elements of 
a liquid helium cryostat with minimum cross-section (for minimizing the heat transfer) and proper 
mechanical resistance. Some examples of numerical results and tables are also presented. 

PACS numbers: 44.10.-K 



I. INTRODUCTION 

Although the heat transfer theory is well known for ordinary conditions [l| and a huge number of studies and 
applications may be found, there is still a lack of experimental data and theoretical solutions for the problem of the 
heat transfer in cryogenic systems. Taking into account the continuous development of the cryogenic techniques Q 
and their applications especially in superconductor systems, it becomes more and more important to minimize the 
heat transfer flux and hence the loss of cryogenic liquids in such installations. 

There are two processes that enable the heat transfer to a liquid helium cryostat built as a Dewar flask :conduction 
and radiation. In figure 1 we present a possible structure of the cryostat, and one may see that the conduction process 
appears due to the sustaining elements of the inner flask. In order to minimize the heat transfer towards such an 
element it has to be made using a thermal insulating material with a cross section as small as possible. Unfortunately 
this reduces the mechanical resistance of the supports, especially at cryogenic temperatures [1| . Various materials have 
been experimented and it seems that the austenitic steels are quite adequate for this purpose, despite their relatively 
high thermal conductivity. Since the tests for an optimal structure are expensive due to the liquid helium evaporation 
that inevitably occurs, it would be desirable to numerically estimate the heat transfer for various materials, shapes 
and dimensions, in order to avoid as much as possible the experimental optimization of the system. 

However, the usual formulae for heat transfer in rods and for their mechanical characteristics fail in cryogenic 
regime, due to the variation with the temperature of the conductance 0| and tensile strength of the materials [|[. 
The aim of this paper is to include as accurate as possible these variations in the formulae necessary for the optimal 
design. 
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FIG. 1: A typical Dewar flask for liquid helium 
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II. HEAT TRANSFER IN A VARIABLE PROFILE CANTILEVER UNDER CRYOGENIC 

CONDITIONS 

We are interested in calculating the rate of heat transfer in a rod placed in vacuum and with the extremities kept 
at two different constant temperatures, T% = 4.2K and X2 = 77K in steady state conditions. 

The heat transfer rate may be described by the Fourier law which states that this is proportional to the gradient 
in temperature and to the normal area through which the heat is flowing 

J = -KSVT (2.1) 

where: 

• J is the rate of heat transfer (heat flux)[Wm~ 2 ] 

•K is the conductivity of the material (constant) [Wm _3 K _1 ] 

• S is the cross section of the rod (constant) [m 2 ] 
•VT is the gradient in temperature [Km -1 ]. 

As there is no conduction heat transfer at the lateral surface of the rod placed in vacuum, the problem is one- 
dimensional and the solution may be written as 

where 

•Ti and T2 are the constant temperatures at the rod's edges [K] 
•L is the length of the rod [m] 

This equation is valid only for constant cross section and conductivity of the rod and small temperature difference. 
If this is not the case, one must divide the whole rod in thin slices to ensure the validity of the eq. (|2.2p . Thus, we 
may write: 

J = K ■ S(x) ■ (2.3) 
dx 

Taking also into account the dependence of the material's conductivity with the temperature: 

K = K(T{x)) (2.4) 

the equation (|2.3I ) becomes 

J = K(T(x))S(x)^ (2.5) 
dx 

For cryogenic temperatures (lower than the Debye temperature), the mean free pass of the phonons may exceed 
the dimensions of the rod and hence the conductivity will be a function only of the specific heat capacity Cy [3 

K{T)= X Cv{T) (2.6) 

where the constant \ is a characteristic of the material. 

According to the classical Debye model [5j , in this regime the specific heat capacity increases with the temperature 
as [| 

C v = aT 3 (2.7) 

involving another material's constant a. 

More general, the Rosseland approximation Q Q allows considering also an additive constant b due to the molecular 
conductivity 



C v = b + aT 3 



(2.8) 
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As for 6 = the Debye dependence is found, we shall consider this more general formula with the parameter b more 
or less significant for various materials. 
From eqs. (l2l)]) . fljjlj) and ([2~8]l we obtain: 



J = X [b + aT 3 {x)}S(x) 



dT{x) 
dx 



Separating the variables we may now integrate this ordinary differential equation 



(2.9) 



r dx 

W) 



X [b + aT 3 (x)]dT(x) 



(2.10) 



(2.11) 



where we took into account that the thermal flux J is constant in steady-state conditions. 
From equation (|2.11|) we obtain the heat transfer rate J as: 



J = x {b(T 2 -T 1 ) + -(Ti-T?) 



dx 



(2.12) 



We may also obtain the temperature distribution along the rod, which will be necessary for mechanical resistance 
calculations. From eq. (|2.5j) we may write 



dT{x) = J- 



dx 



K(T(x)) ■ S(x) 

At a distance I from the T\ source, the temperature T(l) is obtained by integrating eq. (|2.13p 



(2.13) 



T(J) 

/ 

Ti 



x[b + aT 3 (x)]dT(x) = J 



dx 



(2.14) 



X {b[T(l)~T 1 ] + ^[THl)-T^]} = J J 



dx 



(2.15) 



Using eq. (|2.12j) for the heat transfer rate J, the temperature distribution is given by the following equation, that 
may be solved numerically for obtaining the temperature distribution in the rod 



bT(l) + -T\l) = m + -Tt + [b(T 2 - Ti) + i (T 2 4 - Tf. 



dx 



dx 



(2.16) 



III. CONSIDERATIONS CONCERNING THE PROFILE OF THE CANTILEVERS IN CRYOGENIC 

CONDITIONS 



It is known that the general Euler-Bernoulli beam theory neglects the share deformations and the beam deflection 
is given by the equation 
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d 2 d 2 u 

where: 

•E is the you Young modulus 
•J is the second moment area 
•Fis the distributed force (force per length) 

The more general Timoshenko beam theory includes the shear forces leading to a coupled linear partial equations: 
where: 

•p is the density of beam material; 

•Ais the cross-section area; 

•G is the shear modulus; 

•k is the Timoshenko shear coefficient; 

The tensile stress in the beam at the distance y from the neutral axis and in the same plane with the beam and 
the applied force is 

d 2 u My 

° = Ey d^ = T (3 - 4) 
where M is the bending moment proportional with the distance to the applied force _F(figure 2) 

M = F(L - x) (3.5) 

Obviously the tensile stress must not exceed the tensile strength a max which is a characteristic of the material, and 
in cryogenic regime depends on the temperature 
Thus, from eqs. (|3T|) and ([3~5]) it follows 

V -^F{L -x)< a max (3.6) 

One may see that the tensile stress increases linearly with the distance to the applied force, being maximum at the 
fixed end of the cantilever and zero at the mobile end. So, a higher mechanical resistance is necessary towards the 
fixed end, meaning that the cantilever has to be ticker as x decreases. We recall that the heat loss is increasing with 
the cross section of the cantilever, so that it must be kept as low as possible. 

Obviously, the solution is to use a variable cross section area cantilever, thinner in the low tensile effort regions and 
thicker in the high tensile effort ones. The ideal dependence of this cross section should keep the same ratio between 
the allowed tensile stress and the actual one in every point among the beam. This will influence also the thermal 
conductivity of the cantilever and the heat transfer flux [9[. 

For a variable cross section area of the beam, the second area moment area / and y ma x are also variable. 

For finding the optimal profile, we must first choose the shape of the cross-section and find the dependence of the 
second moment area on its geometrical characteristics. 

In most cases, the applied force is in a fixed plane and an I shape cantilever provides the largest second moment 
area with a minimum cross section area. However, in our case if the flask is inclined, the applied force changes the 
relative plane and the I shape cantilever has a much lower admissible tensile stress in such situation. 

That is why a circular or a rectangular shape of the cross section may be a proper choice, despite its lower 
performance in the one plane case. 

a) Circular shape cross section cantilever 

The second momentum area in this case is 

4 

i = -r (3-7) 
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From eq. (1.21) with ?/ max = r we obtain 



and the profile of the cantilever is given by 



AF 

— (L-x) < a max {T(x) (3.8) 



r(x)>[ ^-—{L-x)$ (3.9) 

Of course, one should choose a provision for the parameter r (the radius must have a nonzero value ro at the mobile 
end) so that the recommended equation is 

AF i 

r(x)=r + [ -— -(L-.t)]3 (3.10) 

7ro- max (T(a;)) 

In figure 2 is represented the recommended dependence of the radius on the distance to the fixed end. The most 
difficult part in this equation is generated by the dependence with the temperature of the maximum tensile strength 
of the material. Qualitatively, it is known that it decreases with the temperature decrease in cryogenic regime, but 
an exact formula for this dependence is not available. The problem may be solved only numerically or using some 
approximation. For numerically solving the whole problem, one must use a table containing the values of the tensile 
strength of the material for several temperature within the desired range. Now an approximating polynomial function 
may be obtained with high accuracy, using for example the method of Lagrange interpolation 

71 n TV ^ T 1 

&mBx(T(x)) = ]T C max (Tj) J] AJ_=- j (3.11) 

i=l j=l 1 3 

We plug in this function in eq. p.lOp and together with eq (12.12j) we obtain a coupled nonlinear equations that 
should be solved numerically. 

A simpler but less accurate solution is to approximate the dependence of the tensile stress with a first degree 
polynomial in variable x 



CT m ax(r(a;)) = cr - cx (3-12) 

where the constant c may be obtained as the slope of a line which best approximates the data in the mentioned table. 
Thus, eq. (I3.10|) becomes 

AF i 

r = r + [- AL-x)]* (3.13) 

7T((7o — CX) 

Inserting this equation in (|2.12|) we obtain the following expression for the heat transfer rate in a cantilever with 
circular cross section and optimal profile 

J = X [b{T 2 -T 1 ) + - A {T$-Tt)]h (3.14) 

-l 

has an exact analytical solution that may be expressed in terms of Gauss hypergeometric 
functions 2-F1 (a, 6; c; z) 10] 



where I\ = 



dx 
Si(x) 
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b) Rectangular shape cross section cantilever 

In the case of rectangular cross section, the second moment of area is 

, . wh(x) 3 

m = (3.16) 

where w is the width (horizontal dimension) and h is the height (vertical dimension) of the cantilever. Since the h 
dimension has a greater influence to the overall resistance, only it will depend on x. 
Plugging in this formula in eq. (|3.6|) and considering y max (x) — ^4^, we obtain 

h(x) = [J^L(L-z)]§ +h (3.17) 
wcr max 

Again, taking into account the dependence on the temperature of the strength of the material, we may use the 
Lagrange interpolation formula (|3.1ip for a polynomial approximation, or be satisfied with the linear formula .In the 
last case, the profile equation becomes 



h(x) = ho + 



6F 
w((Jo — bx) 



{L-x)Y 



(3.18) 



Also, inserting this equation in (|2.12[) we obtain the following expression for the heat transfer rate in a cantilever 
with rectangular cross section and optimal profile 



J = X {b(T 2 -T 1 ) + -(T*-T?)]I 2 



(3.19) 



where I? = 



dx 
rh (.(■ ) 



may be numerically calculated. 



IV. NUMERICAL RESULTS AND CONCLUSIONS 



Using our analytical formulae for the cross sections of the cantilever for an optimal mechanical resistance and 
minimal conduction we get the numerical numerical results in Table HI We provided nonzero value at the mobile end, 
where the moment is theoretically null, to prevent the shearing deformation that was not included in the Bernoulli- 
Euler equation. No variation of the tensile strength with the temperature has been considered, so that higher values for 
the radius should be taken towards the mobile end. It suggests that a conical profile would be a good approximation 
for the optimum compromise between mechanical and thermal characteristics of the cantilever. 



The proposed formulae may be used for a proper design of the liquid helium recipient in order to minimize the heat 
transfer and preserve the mechanical characteristics of the system. They include the conduction and tensile strength 
dependance on the temperature in cryogenic regime which may be theoretically estimated or experimentally deter- 
mined for the material of choice. Although they lead to a cumbersome profile of the cantilever, a good approximation 
with a conical one may be considered for practical purpose, assuming that the actual cross section is grater or equal 
to that predicted by our formulae in every point. 
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TABLE I: The longitudinal variation of the radius (in mm) of an austenitic steel cantilever necessary for various weights of the 
liquid helium flask 
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FIG. 2: The optimal cantilever profile for a circular cross section 



[3] G. Claudet , Superfluid helium from physics laboratory to industry, Proceedings of ICEC19, Gistau-Seyfert ed., Narosa 

Publishing House, New Delhi: 743-750 (2002); 
[4] D. G. Cahill and R. O. Pohl, "Lattice Vibrations and Heat Transport in Crystals and Glasses," Ann. Rev. Phys. Chem. 

39: 93-121, (1988); 

[5] P. Debye,Zur Theorie der spezifischen Waerme, Annalen der Physik (Leipzig) 39(4), p. 789 (1912); 
[6] S. Rosseland, Ap Phys J, 61, 424,(1925); 

[7] S. Rosseland, " The Principles of Quantum Theory" , Ergnzungsband, pp. 243-249, (1936); 

[8] M. R. Sridhar and M. M. Yovanovich, "Thermal Contact Conductance of Tool Steel and Comparison with Model," Int. J. 

of Heat Mass Transfer 39 (4): 831-839, (1996); 
[9] M. M. Yovanovich, P. Teertstra, and J. R. Culham, "Modeling Transient Conduction From Isothermal Convex Bodies of 

Arbitrary Shape," Journal of Thermophysics and Heat Transfer 9 (3): 385-390, 1995; 
[10] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, 1965. 



